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A Note on the landing of a 
Cracked Strip 

by 

Rag ip Erdol and F. Erdogan 
, Lehigh University* Bethlehem, PA 

Abstract. The objective of this note is to describe a technique for 
calculating the stress intensity factors in a strip under bending by 
treating the strip as a two-dimensional continuum rather than a simple 
beam in evaluating the crack surface tractions used for the, solution 
of the perturbation problem. 

1 . INTRODUCTION 

A long strip or beam containing a through crack perpendicular to 
the sides has been one of the most widely studied problems in linear 
fracture mechanics (see, for example, references [1-9]) • The import- 
ance of the problem lies in the fact that its geometry approximates a 
very common structural component and a standard test specimen. In the 
case of three or four point bend tests the specimen has an edge crack. 
Generally, in the existing solutions involving bending it is assumed 
that the uncracked beam is under a linear stress distribution, thus 
ignoring the perturbations caused by the supports and the loading 
fixtures. 

The main objective of this note is in a "long" beam to examine 

the deviation from linearity in the stress field for a given state of 

^*^This work was supported by NASA-Langley under the Grant NGR 39-007- 
011 and by NSF under the Grant ENG 73-045053 AOl. 
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external loads and to study the effect of this deviation on the stress 
intensity factors for internal and edge cracks. 

2. STRESSES IN UNCRACKED BEAM 

Consider the infinite strip or the beam show in Figure 1. Let 
the external forces be 


-Oyy(x,0) = p^(x) = p^(-x) , a*y(x,0) = 0 , (1) 

-o^y{x,h) = pg(x) = P2(-x) , -cjj^y{x,h) » 0 . (2) 


where h is the height of the beam p-| and P 2 are known functions 
satisfying 


/OO fOO 

J P^(x)dx = J p 2 (x)dx 


(3) 


Ignoring the crack and using the standard Fourier transform technique 
(see, for example, [10]), after some routine manipulations the stress 
component of primary interest may be expressed as 


9 r 

= tF J Sj^(a,y)cosax da 


(4) 


where 

i>x(a»y) = {(l-2a^h^-2ah-e^“^)S^(a) 

+ [0+3ah)e“^ - (l+ah)e"“^]S 2 (a) 
+ ay(-l+2ah+e^“^)S^(a) 
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- ay[(l+2ah)e“^ - 

+ ^ {(l+2ah-2a2h2 - e"2“^)S^{a) 
+ (:(-l+ah)e“^ + (l-3ah)e"“'^]S2(a) 

+ ay(l+2ah-e'^“*^)S^(a) 


+ ay[»e“^ + (l-2ah)e’“^]S2(a)} , 

(5) 

D(a) = - 4a=^h^ - 2 , 

(6) 

S^(a) “ j^p^(x)cosax dx , 

(7) 

S 2 (a) “ j P 2 (x)coxax dx . 

(8) 

3. THE CRACK PROBLEM 


To obtain the solution of the crack problem shown in Figure 1 

under the external loads given by (1) and (2) one has to superimpose 

on the solution obtained in the previous section for the uncracked 
a 

strip^disturbed stress state found from a cracked strip in which the 
following self-equilibrating crack surface tractions are the only ex- 


ternal loads: 

. <^xy^0»y) = 0 , a<y<b , (9) 
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where oj is given by (4). The formulation of this problem too is 
relatively straightforward. Following, for example, the technique 
described in [6] or [7], the problKii may be reduced to the following 
integral equation: 

fa ►^(y»t)]f(t)dt = -TT^oyo.y) , a<y<b , (10) 

subject to 

Lf(t)dt = 0 (11) 

where the input function is obtained from (4), y is the shear 
modulus, k=3-4v for plane strain (e.g., cylindrical bending of a 
plate with a long surface crack), k=(3-v)/(1+v) for plane stress 
(e.g., a beam), v is the Poisson's ratio, crack extends along the y 
axis from a to b, the unknown function f is defined by 

f(y) = ^ U(0.y) . (12) 

u is the X-component of the displacement vector, and the kernel k(y,t) 
is given by 

k(y.t) = r [(e“‘^-e““^)"-4a2h"]“'' {[-l-aab+e^"*^ 

Jo 

-2at(e2“''-l)][| + c(h-y) 

+ [2ah-He"^“'’-4tha2][- |+ a(h-y) 
t [l+2ah-e"*“*'+2ci(h-t)(e^“''-l)][|+oy 


> 4 . 


( 13 ) 


+ 1 t [2ah-1*e' 


2ah 


- 4(h-t)a^h][-2- ay + ^ }c|t>< , 

For an internal crack 0<a<b<h, the kernel k(y,t) is bounded in 
the closed domain a<(y,t)<b, and the integral equation may be solved 
quite simply by using the technique described in [11]. The stress in- 
tensity factors are then defined and calculated from 


k{a) = lim /2(a-y) S » 

k(b) = lim /2Ty-b) a^^(y.O) = -lim /2(b-y) ^ f(y) . 


y-^-b 


y^b 


(U) 


(15) 


However, in the case of an edge crack 0=a<b<h, the kernel k(x,t) 
is no longer bounded in the closed-domain 0<(y,t)<b. The singular 
part of k which becomes unbounded as y and t approach the end point 
a=0 may be separated by considering the asymptotic behavior of the 
integrand in (13). Expressing 


k(y.t) = kg(y,t) + kf(y,t) 


(16) 


where k^ is bounded in a<(y,t)<b, the singular part kg is found to be 

k_(y»t) = [ (2-ay-3at+2yta^)e"^^^^'^^da 
^ Jo 


1 




4y^ 


t+y Tt+y7^ ' (t+y) '‘ ‘ 


(17) 


Noting that the integral equation (10) is still valid for a=0, after 


- 5 - 


I I - 1 I .. .. 1 / \ 


separating it may be solved by following the procedure described 
In [63. 

4. AN EXAMPLE AND SOME RESULTS 

As an example consider the beam supported at two points x=+d on 
the side y=0 and subjected to uniform pressure on y=h, -c<x<c (see 
Figure 1). In reality all loading problems involving solid grips and 
fixtures are contact problems. For example, the supports at x=+d, y=0 
are generally rollers often with different elastic properties. In the 
three or four point bending problems the same is true for the loading 
fixtures on y=h. In this case the contact areas and the distribution 
of contact pressures as well as the crack surface displacement are un- 
known and the problem is a coupled contact-crack problem which is some- 
what more difficult to formulate. However, if the contact areas are 
small and are sufficiently far from the crack region, one may replace 
the unknown contact stresses with statically equivalent known trac- 
tions gaining considerable simplicity in the formulation of the problem 
without sacrificing too much in the accuracy of the results. Thus, in 
the present example it will be assumed that the reactions at the sup- 
ports~y=0, x=+d are concentrated forces of magnitude P (per unit thick- 
ness). In the three point bending problem the pressure under the loading 
pin will be assumed to be constant with c=0.05h. The external loads 
and their Fourier transforms may then be expressed as follows (see 
equations 1 to 4): 

p^(x) = P6(x-d) , (18) 
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a 

[p/c , 0<x<c , 

(19) 

U 

Pp(x)»< 

1 0 » C<x«» , 


S-|(a) = P cos (da) , 

(20) 


r P sin(ca) .. 

- c a - ' ' 

(21) 


In evaluating the input function o' from (4)-(8) one runs 

A A 

into 


convergence difficulty for values of y around 0 and h. To avoid this 
difficulty the damping of the integrand as a->~ is increased by separ- 
ating the asymptotic value of S {a,y) and integrating it in closed 

A 

form. Equation (4) may then be expressed as 




f” 

[S„(a,y) - S^(a,y)]cosax da + sjx.y) , 

^0 

(22) 


S„(a,y) = (ay-l)e"'^^S^(a) - (ay+l-ah)e"“^^'^^S 2 (a) , 

(23) 

‘-tf 

=■ 1 

r 

1 s (a,y)cosaX da , 

(24) 


s^(0,y) = P 

y ^ v(v^-d'') ^ h-y 1 t c ^ 

• 



.. — 

(25) 


The stress distribution o' (0,y) obtained from the elasticity 
theory as described by equations (20) to (25) and that obtained from 
the simple beam theory for the loading condition shown in Figure 1 is 
given in Table 1. The normalization stress o^ used in this and in the 
subsequent tables is the surface stress given by the beam theory which 
may be expressed as 
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°o ~ ^ • (26) 

The table shows that in the elasticity results there is some 
deviation from that of linear beam theory. The indication is that 
the sign as well as the magnitude of the deviations depend on details 
of the loading condition. However, at least for the particular load- 
ing considered here, the relative magnitude of the deviation does not 
appear to be very significant (approximately 6 percent at y=0). 


Table 1. The stress profile c;iY(0, y)/aQ for the uncracked beam 
calculated from (22-25) for various dimensions c and d. 
The normalization factor: o = 6P(d-c/2)h^. 


y/h 

c/h = l 
d/h=l 

c/h = 0.05 
d/h =2 

c/h = 0.05 
d/h =4 

Beam Theory 
(h/2-y) 

0 

1.0589 

0.9435 

0.9399 

1 

0.1 

0.7949 

0.7276 

0.7384 

0.8 

0.2 

0.5595 

0.5363 

0.5491 

C.6 

0.3 

0.3486 

0.3615 

0.3680 

0,4 

0.4 

0.1566 

0.1971 

0.1921 

0.2 

0.5 

-0.0233 

0.0386 

-0.0192 

0 

0.6 

-0.1993 

-0.1179 

-0.1528 

-0.2 

0.7 

-0.3801 

-0.2767 

-0.3259 

-0.4 

0.8 

-0.5750 

-0.4452 

-0.5039 

-0.6 

0.9 

-0.7938 

-0.6723 

-0.7109 

-0.8 

0.95 


-.09153 

-0.9991 



-0.9204 



-0.9 


Table 2 gives the stress intensity factors for the internal crack 
shown in Figure 1. In this example the crack tip a was fixed at a=0.1h 
and b was varied until the stress intensity factor at b, k(b) became 
(approximately) zero or negative. Note that for a very small crack 
(i.e., a=0.1 h, b=0.1001h) k(a) 2 k(b) = a^^(0.0.1h)/r which is the 
expected result for an infinite plane with line crack of length 2)l=b-a 
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Table 2. Stress intensity factors in a beam with an internal 
crack under the loading condition shown in Figure 1 
cf = 6P(d-c/2Vn^, £= (b-a)/2, and a=0.1h=constant. 
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pressurized by (see the rows y/h=0.1 in Table 1 and b/h=0.1001 in 
Table 2). For fixed a/h the stress intensity ratio k(a)/aQ/4 remains 
approximately constant, meaning that for increasing crack length, k(a) 
itself increases as /i. Also note that k(b) is reduced rapidly as the 

crack tip b moves towards and into the compression region. For the 

beam theory extensive results are given in [7]. Table 2 shows the re- 
sults from [7] which coincides with the parameters a/h and b/h in the 
table. 

In the important case of the edge crack (i.e., a=0, b<h in Fig- 
ure 1) the results for c=C).05h and d=4h are shown in Table 3. The 
stress intensity factor obtained from the beam theory is reproduced 
from [7]. For a very short crack one would expect to recover the 
result of the semi -infinite plane with an edge crack of length b pres- 
surized by 0 * (0,0) given in Table 1. Indeed for b/h=0.001 it is seen 
that 

k(b) _ 1.0537 . , ..o.,, 

ci^(0,0)y^ “ 0.9399 " 

AX 

which is the half plane result/*^ 

The results shown in Tables 2 and 3 indicate that, depending on 
the loading condition, the elasticity solution may be different than 

^*^For higher values of b/h the convergence of the numerical analysis 
giving the stress intensity factor is somewhat slow. The results 
shown in the table are obtained by letting k(b)“A/n^+B/n+C and 
extrapolating the results to n=«, where n is the number of collo- 
cation points in the solution of the related integral equation. 
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Table 3. Stress intensity factor in a beam with an edge crack under 

the loading condition shown in Figure 1 {c=0.05h, a=0, b<h,c|R^h), 
a^=6P(d-c/2)/h^ 


b/h 

k(b)/ao/P 

Elasticity 

Beam Th. 

0.001 

1.0537 

1.12 

0.1 

0.9731 

1.05 

O.F 

0.9750 

1.06 

0.3 

1.0375 

1.12 

0.4 

1.1661 

1.26 

0.5 

1.3907 

1.50 

0.6 

1.7839 


0.7 

2.5127 



the beam solution. As to how significant this difference is depends on 
the degree of accuracy required in the particular application as well 
as the details of the loading condition. It should only be pointed out 
that the procedure outlined in this paper is very straightforward and 
would give the results for any given state of loading to any desired 
degree of accuracy. 


P 
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